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Abstract. Class objects, which are thought to be the youngest protostars, are identified in terms of NIR or radio emission 
and/or the presence of molecular outflows. We present combined hydrodynamic and radiative transfer simulations of the 
collapse of a star-forming molecular core, which suggest two criteria for identifying dense cores with deeply embedded very 
young protostars that may not be observable in the NIR or radio with current telescopes. We find that cores with protostars 
are relatively warm (T> 15 K) with their SEDs peaking at wavelengths < ITO/xm, and they tend to appear circular. 
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1. Introduction 

Class objects, the youngest protostars identified so far, are 
distinguished from prestellar cores, their precursors, based 
on 2 criteria (see Andre et al. 2000), (i) the presence of a 
NIR source in the centre of the core, or (ii) the presence of 
compact radio emission or molecular outflows. These criteria 
depend on the sensitivity of the telescope used. Very young 
protostars are difficult to observe in the NIR, because they 
are deeply embedded. Additionally, the amount of their radio 
emission is uncertain (see Gibb 1999), and probably very low 
(Neufeld & Hollenbach 1996). Thus, a null detection with 
a specific telescope does not necessarily mean that a source 
does not exist. This was demonstrated recently with the core 
L1014 which was thought to be prestellar because no source 
was detected by IRAS, but later was suggested to be proto- 
stellar after a NIR source was detected by the more sensitive 
SPITZER (Young et al. 2004). Hence it is important to con- 
sider ways in which the presence of young embedded pro- 
tostars might be inferred indirectly. In this paper we present 
SPH simulations of the collapse of a turbulent molecular core, 
combined with 3 -dimensional Monte Carlo radiative trans- 
fer simulations at different stages of this collapse. The sim- 
ulations predict dust temperature fields, SEDs, and isophotal 
maps of very young protostars that are deeply embedded in 
their parental cores, and provide criteria for distinguishing 
between genuine prestellar cores and cores that contain very 
young protostars (Stamatellos et al. 2005). 



2. The model 

We use the SPH code Dragon (Goodwin et al. 2004) to sim- 
ulate the dynamics of a collapsing molecular core and the 
Monte Carlo radiative transfer code Phaethon (Stamatel- 
los & Whitworth 2003, 2005; Stamatellos et al. 2004) to treat 
the radiative transfer within the core, at different time-frames 
during its collapse. 

2.1. The collapse of a molecular cloud 

The initial conditions in the core, before the collapse, are set 
according to observations (e.g. Kirk et al. 2005). Prestellar 
cores have typical extents from 2000 to > 15000 AU, and 
masses from 0.05 to 10 M©. Their density profiles are flat in 
their centre and fall off as r~^, where n = 2 — 4. These fea- 
tures are well represented by a Plummer-like density profile 
(Plummer 1915), 



(l + (r/rJ2)2 



(1) 



and falls 



The density profile is almost flat for r < 
off as in the outer envelope (r > r^). We set = 



3 X 10-^^ 



g cm 



5,000AU, andr^ = 50,000AU. 
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The core mass is then 5.4 Mq. We impose a low level of 
turbulence (a^u^j3 = ^turb/I^grav I = 0-05). 

The collapse of the turbulent core leads to the formation 
of a single star surrounded by a disc, and the material in the 
disc then slowly accretes onto the star. A sink particle, rep- 
resenting the star and the inner part of the disc, is created 
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3. 1 Dust temperatures 



Table 1. Model parameters 



id 


time (yr) 


M.(M0) 


Rs (AU) 


Zs (AU) 


Ms(M0) 


M^(M0/yr) 


Ltot(L0) 


T.(K) 


description 


tl 


5.3 X 10^ 
















prestellar core 


t3 


6.0 X 10^ 


0.20 


4.0 


0.4 


0.09 


1 X 10"^ 


27.2 


7650 


Class protostar 


t5 


6.9 X 10^ 


0.53 


4.0 


0.4 


0.01 


4 X 10"^ 


2.5 


4200 


protostar is off centre 



Zs, Rs- Smartie dimensions M^: Stellar mass T^: Star temperature 

Ms: Smartie mass M^: Accretion rate onto the central star Ltot: Total star luminosity (intrinsic +accretion) 



where the density first exceeds Pcmx = 10~^^ gcm~^. We 
use a newly developed type of sink called a smartie, which is 
a rotating oblate spheroid. The star is a point mass at the cen- 
tre of the smartie.The smartie gains mass from the accreting 
material and can lose mass by ejecting material in the form 
of protostellar jets. The jets are simulated by SPH particles 
which are created at rate O.lTM^/t^^^^ and launched along 
the rotation axis with speed lOOkms"^ to form a bipolar 
outflow. These jets push the surrounding material aside and 
create hourglass cavities about the smartie rotation axis (see 
Stamatellos et al. 2005 for details). 

2.2. Radiative transfer 

We perform radiative transfer simulations on 3 time-frames 
during the collapse of a star-forming core (see Table [ij. We 
focus our attention just before and just after the formation of 
the first protostar in the core. 

For the radiative transfer simulations, we use a version of 
Phaethon, a Monte Carlo radiative transfer code which we 
have developed and optimised for radiative transfer simula- 
tions on SPH density fields. The main new features of the 
optimised version are (i) that it uses the tree structure in- 
herent within the SPH code to construct a grid of cubic ra- 
diative transfer (RT) cells, spanning the entire computational 
domain, the global grid; and (ii) that in addition it constructs 
a local grid of RT cells around each star, a star grid, repre- 
senting a flared disc, to capture the steep temperature gradi- 
ents that are expected in the vicinity of a star (Stamatellos & 
Whitworth 2005). 

The core is illuminated externally by the interstellar radi- 
ation field and internally by the newly formed protostar (once 
it has formed). For the external radiation field we adopt a re- 
vised version of the Black (1994) interstellar radiation field, 
which consists of radiation from giant stars and dwarfs, ther- 
mal emission from dust grains, cosmic background radiation, 
and mid-infrared emission from transiently heated small PAH 
grains (Andre et al. 2003). The protostar luminosity is domi- 
nated by the luminosity due to accretion since the mass of the 
young protostar is small (< 0.6 Mq) and the accretion rate is 
highO lO-^M0/yr). 

3. Characteristics of prestellar cores and young 
protostars 

3.1. Dust temperatures 

In Fig. [I] we present the density fields and dust- temperature 
fields for the 3 time-frames in Table [l] We plot the density 




Fig. 1. Cross sections of density and dust temperature on a 
plane parallel to the x = plane of 3 different time-frames 
(first row: collapsing prestellar core; second row: Class 
object; third row: Class object). 



and the temperature on a cross section through the computa- 
tional domain parallel to the x = 0. This plane passes through 
the protostar, or, if there is no protostar (as in the first time- 
frame), through the densest part of the core. The plots show 
only the central region of the core (5000 AU x 5000 AU). 

Our results for the temperature are broadly similar to 
those of previous ID and 2D studies of prestellar cores and 
Class objects (Evans et al. 2001; Zucconi et al. 2001; 
Shirley et al. 2002; Stamatellos & Whitworth 2003). The 
prestellar core is quite cold (5 to 20 K). As soon as a pro- 
tostar forms, the region around it becomes very hot (up to 
the dust destruction temperature), but the temperature drops 
below ~ 20 K beyond a few 1000 AU from the protostar, be- 
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Fig. 2. SEDs of 3 different time-frames from 6 different angles (the SED is weakly dependent on viewing angle so that some 
of the curves overlap; the flux is higher when the protostar is viewed pole-on, i.e. through the hourglass cavity that is created 
by the jets), (a) first column - collapsing prestellar core. The SED peaks at ~ 190 /im (temperature T^^^ ~ 13 K). (b) second 
column - Class object. The emission of the system peaks at A ~ 80 to 100 /im (depending on the observer's viewing angle) 
(^EFF ~ 25 to 31 K). (c) third column - Class object. The protostar has moved out of the central region by ~ 100 AU, and 
the attenuated stellar emission can be seen at short wavelengths (1 to 50 /im).The peak of the cloud emission is at A ~ 150 /im 
(T,,, ^17K). 




Fig. 3. 850 /im isophotal maps of 3 different time-frames (first column: collapsing prestellar core', second column: Class 
object', third column: Class object). Cores with protostars (second and third columns) are more centrally condensed than 
prestellar cores (first column). The axes (x, y) refer to the plane of sky as seen by the observer 
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Fig. 4. As in Fig.[3l but after convolving with a Gaussian beam (FWHM = 2300 AU, e.g. 15'' for a core at 150 pc). Cores 
with protostars appear more circular, but otherwise similar to the prestellar core. 



cause of the high optical depth in the dense accretion flow 
onto the protostar. 

3.2. SEDs 

The SEDs of the 3 time-frames in Table [l] are presented in 
Fig. 121 These SEDs have been calculated assuming that the 
core is at 140 pc. SEDs are plotted for 6 different viewing 
angles, i.e. 3 polar angles (^ = 0°, 45°, 90°) and 2 azimuthal 
angles (^ = 0°, 90°). 



The effective temperature of the core, as inferred from 
the peak of the SED, rises and falls with the accretion lu- 
minosity of the protostar. For a prestellar core, the SED 
peaks at ~ 190 /im, implying T^^^ ~ 13 K. Once the pro- 
tostar has formed and inputs energy into the system, the peak 
moves steadily to shorter wavelengths, reaching ~ 80 /im 
(T^pp ~ 31 K) as the accretion luminosity reaches its maxi- 
mum, and then moving back to longer wavelengths again as 
the accretion luminosity declines. By the final frame it has 
reached - 150 /im (T^^^ - 17 K). 
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The peak of the SED of a prestellar core is independent 
of viewing angle, since the core is optically thin to the radi- 
ation it emits. In contrast, the peak of the SED of a Class 
object does depend on viewing angle, albeit weakly, because 
of the presence of an optically thick disc around the protostar. 
However, even allowing for variations in the viewing angle, 
the SED of a Class object does not peak at the wavelengths 
characteristic of prestellar cores (~ 190 /im). 

Our model predicts that a young protostar embedded in 
a core is not observable in the NIR, unless it is displaced 
from the central high-density region. This refers to the very 
early stages of protostar formation. This result contradicts 
the recent results of Whitney et al. (2003) and Young et al. 
(2004); they predict that NIR radiation from a deeply em- 
bedded young protostar is observable. We attribute the dif- 
ferences to the fact that these models use different densities 
for the region near the protostar and different opacities. Ul- 
timately, the density distribution within 500 AU of the the 
centre of a core or a protostar is not well constrained, and 
therefore radiative transfer calculations have to rely on theo- 
retical models of core collapse. Dust opacities are also poorly 
constrained (e.g. Bianchi et al. 2003). 

4. Isophotal maps 

In Fig.[3lwe present isophotal maps at 850 /im for the 3 time- 
frames in Table 1. At 850 /im the core is optically thin, and 
the temperature does not vary much (~ 10 to20K, apart 
from the region very close to the protostar in time-frames t3 
and t5), so the maps are effectively column density maps. 
Class objects are more centrally condensed than prestel- 
lar cores. They also show more structure, due to bipolar out- 
flows, which clear low-density cavities. 

In Fig.|4|we present the isophotal maps after convolving 
them with a Gaussian beam having FWHM = 2, 300 AU. 
Assuming the cloud is 150 pc away, this corresponds to an 
angular resolution of 15'^ which is similar to the beam size 
of current submm and mm telescopes (e.g. SCUBA, IRAM). 
On these maps. Class objects look very similar to prestel- 
lar cores, apart from the fact that they tend to appear more 
circular and featureless. 

The reason for this is that the emission from a core that 
contains a protostar is concentrated in the central few hundred 
AU, and so, when it is convolved with a 2,300AU beam, 
it produces an image rather like the beam, i.e. round and 
smooth. In contrast, the emission from a prestellar cores has 
structure on scales of several thousand AU, much of which 
survives convolution with a 2, 300 AU beam. Thus we should 
expect cores that contain protostars to appear rounder than 
prestellar cores, and indeed this is what is observed (Good- 
win et al. 2002). 

5. Conclusions 

Our simulations indicate that a newly-formed protostar may 
not be observable in the NIR with current telescopes, because 



it is too deeply embedded. This result reflects the high den- 
sities which our model predicts in the immediate surround- 
ings of a newly-formed protostar. Thus, it is suggested that 
cores with no observed NIR radiation may in fact contain the 
youngest protostars. 

Based on the results of the combined hydrodynamic and 
radiative transfer simulations, we propose two criteria for 
identifying cores which - despite appearing to be prestellar 
- may harbour very young protostars: 

(a) They are warm (T > 15 K) as indicated by the peak 
of the SED of the core (Apeak < 170 /im). This criterion re- 
quires that the peak of the SED can be measured to an accu- 
racy of ~ 30 /im, which should be feasible with observations 
in the far-IR from ISO and the upcoming Herschel mission, 
and in the mm/submm region from SCUBA and IRAM. We 
have presumed that the core is heated by the average interstel- 
lar radiation field, and the criterion would have to be modified 
for a core irradiated by a stronger or weaker radiation field. 

(b) Their submm/mm isophotal maps are circular, at least 
in the central 2000 to 4000 AU. 

These criteria can be used to identify "warm", circular 
cores that may potentially contain the youngest protostars. 
These cores can then be probed by deep radio observations 
with the VLA. Recently, we carried out 3.6 cm VLA obser- 
vations for 4 of these cores and results will be presented in a 
future publication. 
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